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It is important to study the behavior of a i-error correcting quantum 

code when the number of errors is greater than t, because it is likely that 

there are also small errors besides t large correctable errors. We estimate 

the fidelity of a i-error correcting stabilizer code over a general memoryless 

channel, allowing more than t errors. We also show that the fidelity can 

be made arbitrary close to 1 by increasing the code length. 

PACS numbers: 03.67.Hk, 89.70,+c 



I. INTRODUCTION 

In the study of the quantum error-correcting codes, it is usually assumed that only a 
small number of qubits are affected and the rest of qubits are left unchanged. However, it 
is important to study the behavior of a t-error correcting quantum code when the number 
of errors is greater than t, because it is likely that there are also small errors besides t large 
correctable errors. The goal of this paper is to provide a lower bound for the fidelity of the 
quantum error correction under the general noise model without any approximation. The 
fact that quantum error-correcting codes work under the general noise model seems a folklore 
result, and the original contributions of this paper are a rigorous proof and a quantitative 
relation between the fidelity and the noisiness of the channel. 

The following researches have been done prior to this paper. It has been informally argued 
in Ref. Sec. VI] that those small errors do not result in a large error in the recovered 
quantum state. The first rigorous analysis was done in Ref. P, Sec. 5.4], in which the 
authors assumed that the channel was memoryless, that is, each qubit interacts with different 
environment, and there was a scalar multiple of the identity operator in an operator sum 
representation of the channel superoperator. In Ref. || Sec. 7.4.2], quantitative relations 
between the fidelity and the noisiness were given for two specific classes of memoryless 
channels. Aharonov and Ben-Or |4], Sec. 8] analyzed the fault-tolerant quantum computation 
under the general noise model that is equivalent to a memoryless channel, and showed that 
if the channel is not too noisy then the error-free computation is possible. However, they did 
not provide a quantitative condition of general channels allowing the error-free computation. 
(They provided that of restricted channels.) 

In this paper we assume that a unitary representation of the channel superoperator has 
large identity component (Assumption 3), and we give a lower bound (Eq. (|7|)) for the 
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average of the fidelity between the original state and the recovered state without using 
any approximation, where the average is taken over the measurement outcome in the error 
correction process. As a consequence we show that the average of the fidelity can be made 
arbitrary close to 1 over a general memoryless channel by increasing the code length. This 
fact has been proved only over specific classes of quantum channels [0, |3| . Our estimation is 
restricted to the stabilizer quantum codes introduced in Refs. || |6], [/]], which include almost 
all good quantum codes discovered so far. It should be noted that the essential idea in our 
analysis already appeared in Ref. || Sec. 7.4]. 

This paper is organized as follows: In Sec. |TJ we introduce notations used in this paper, 
and review the stabilizer quantum codes and their error correction process. In Sec. [Hi] we 



give a lower bound for the fidelity. In Sec. |V| several consequences and generalizations are 
discussed. 



II. PRELIMINARIES 



A. Notations 

Let H be a Hilbert space. We denote by S(H) the set of density operators on H. For 
a density operator p on H and a state vector \ip) G H, the fidelity [§, || between them is 
defined by 

F(hA),p) = <VW>. 

It measures how close \ip) and p are. 

In this paper we consider t-error correcting [[n, k]] binary quantum codes unless otherwise 
stated. Let H2 be the Hilbert space of dimension 2. Let T be a superoperator on H2, that is, 
a trace-preserving completely positive linear map from S(H2) to S(H2). We assume that the 
channel is represented by T, which means that when we send a density operator p G S(H 2 ) 
through the channel we get G S(H 2 ) at the receiving end. The channel considered 
in this paper is assumed to be memoryless. So when we send a state p G ") we get 

r® n ( P ). 

We shall review the unitary representation of a superoperator [[H]]. A simplified proof can 
be found in Ref. Appendix]. Let T be a superoperator on a Hilbert space H. Then there 
exist a Hilbert space He, a state vector \0e) G He, and a unitary operator U on H ® He 
such that 

T(p) = Tt e (U(p®\0e)(0e\)U*), (1) 

for all p G <S(H), where Tr^ is the partial trace over He- That is called a unitary represen- 
tation of r. 



B. Stabilizer quantum codes 



In this subsection we review the method of quantum error correction proposed in Refs. 
& ill- Let 
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and E = {±wi® ■ ■ ■ ®w n }, where Wi is either i, a x , a z or a x a z . The set E is a noncom- 
mutative group with matrix multiplication as its group operation. Let S be a commutative 
subgroup of E. A quantum error-correcting code Q C H®™ is defined as an eigenspace of S. 

For M G E we define MQ = {M\<p) : \f) G Q}- The set MQ is also an eigenspace of S 
for any M G E. Moreover {MQ : M G -E} is equal to the set of eigenspaces of S. It follows 
that every eigenspace of S has the same dimension. Let dimQ = 2 k . Then there are 2 n ~ k 
eigenspaces of S. Let S' = {N G E : MN = NM for all M e S}. It is known that 

S' = {M e E : MQ = Q}. (2) 

We shall describe the error correction procedure. Let H env be the Hilbert space repre- 
senting the environment around the channel. Suppose that we send a pure state \<p) G Q, 
and the environment is initially in a pure state |0 en v) G H env . Suppose also that we receive 
an entangled state \ip) G Hf n <8> H env . We measure an observable of if® n whose eigenspaces 
are the same as those of S. Then the state \ip) is projected to \ip r ) G Q' ® -ffenv, where Q' is 
some eigenspace of S. 

We will define the weight of an operator M G E for error correction. Let M = ±w\ (g) 
■ • where Wi is either /, cr x , a z or a x cr^. The weight of M is defined to be : it), 7^ /}, 

and denoted by We define the numbers d and d' by 

d = mm{w(M) : M G 5' and ± M <£ S}, 
d' = mm{w(M) : M G S' and ± M ^ /}. 

The number d is called the minimum distance of Q. The code Q is s& id to be pure if d = d' 
and impure if d > d' . We define t = [(d — l)/2j . 

There are many operators M e E such that MQ = Q'. Let M be an operator whose 
weight is minimum among them. Note that if the weight of M is greater than \_(d' — l)/2j 
then there may be another operator M' such that w(M') = w(M), M'Q = Q' and M 7^ ±M'. 
We guess that the original pure state is (M -1 ® J env )|^'), where J env is the identity operator 
on H env . 

If the number of errors < t, then is the tensor product of \ip') G Q' and some pure 
state in H env , and M -1 !^') = \(p). However, we do not make such assumption, and we shall 
analyze the closeness (fidelity) between \<p) and M~ 1 Tr^ cnv (|'0 / )('0 / |)(M" 1 )*. 

We shall use the following fact later. 
Proposition 1 Let M' G E be an operator such that M'Q = MQ. If ±M' MS then 
w(M') >t, where MS = {MN : N e S}. 

Proof. If w(M) > t then w(M') > t by the definition of M. Suppose that w(M) < t and 
w(M') < t. Then w(M~ 1 M') < 2t < d, M ■ M' G S' by Eq. (|), and M _1 M' ^ S. This 
contradicts to the definition of d. ■ 



III. LOWER BOUND FOR THE FIDELITY 

In this section we consider the fidelity between the original state and recovered state. Let 
T be the channel superoperator of if 2 as in Sec. [II A| . Since I, a x , a z and a x a z form a basis 



of linear operators on if 2, in a unitary representation of T we can write U in Eq. ([!]) as 



i <g> L ,0 + cr x (g) Li, + (J* <g) L ,i + (JxCr^i,!, 
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where L it j is a linear operator on a Hilbert space H E . Let \0e) be the initial state of H E . 
Lemma 2 We retain notations as above. ||L 0i oO_e|| < 1, where || • || denotes the norm of a 
vector ■. 

Proof. Let {|0), |1)} be the orthonormal basis such that a x |0) = |1), a x \l) = |0), <r 2 |0) = |0), 
and cr 2 |l) = — 11). Then we have 

U(\0) ® |0 S » = |0)®(Lo,o|0 B )+Lo,i|0 s )) + |l)®(L 1 ,o|0 s )+ J L 1 ,i|0 s ), 
C/(|1)®|0 E }) = |l}®( J Lo,Q|0 E )-L 0)1 |0 £ }) + |0)(g)(L 1) o|0 £ }-L 1)1 |0 E ). 

Since both vectors are of unit norm, it follows that 

||^o,o|0b)+^o,i|0b>|| < 1, 
||£o,o|0js) -io,i|0^)|| < 1. 

We conclude ||Lo,oO_b|| < 1- ■ 
Assumption 3 Assume that 

||L ,iO £ || 2 + H^oOeII 2 + H^i^o^ll 2 = p. 

Hereafter we denote H E n by H env and \0e) <£> • • • ® \0e) by |0 env ). Suppose that we send 
\<p) G Q and the recovered state is M- lr Tr Hcnv (\xp') (^^(M" 1 )* as in Sec. [QB[ 

We shall consider the average of .F(|</?), M _1 Tr^ env (|'0 / )('0 / |)(M _1 )*) for an arbitrary fixed 
state G Q under the assumption that the channel is memoryless. The superoperator of 
the channel is r® n . Let Z2 = {0, 1} with the addition and the multiplication taken modulo 
2. For a = (a±, . . . , a n ) G Z?J , we define 

X(a) = 

Z(a) = o-f <g> • • • (8) a a z n . 
Then a unitary representation of T® n can be written as 

£ X(a)Z(6)®L^, 



where 



Let |^) be as in Sec. [II B| . By notations defined so far, can be written as 



|V>> = E X(a)Z(6)|^)®L S6 H0 env ). 

Let Q' be an eigenspace of S. We shall consider the probability Pq> of |^) being pro- 
jected to Q' ® H env after the measurement. Let (agr, 6q/) be a pair of vectors such that 
X{a Q t)Z(b Q ,)Q = Q' and that if M'Q = Q' then w(M') > w(X(a Q ,)Z(b Q ,)). Observe that if 
w(X(aQi)Z(bQi)) > l(d'— l)/2j then there may be another operator M' such that M'Q = Q', 
w(M') = w(X(aQi)Z(bQi)) and M' 7^ ±X (Sqi) Z (bQi) . This implies that (ciq', 6q/) is not 
uniquely determined by Q'. One may choose whichever (a,Q>, 6q/) provided that X{a,Qi)Z(bQi) 
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has the minimum weight (see also Sec. gVg) . Let T Q , = {(c, d) G Z£xZ£ : X(c)Z(cf)Q = Q' 
and ±X(c)Z(d) (£ X(a Q ,)Z(b Q ,)S}. We define 



Wq>) 



E 



±X(c)Z(d)eX(a Q ,)Z(b Ql )S 
(S Q i ,dQi)eT Q i 



(3) 



Observe that 

Q' is an eigenspace of 5 

and the projection of \ip) to Q' ® if e nv is |ctq/ + cTq,). Thus Pq> is given by ||ctq' + cTq,|| 2 . 

Let = |ctq/ + cTq,)/||ctq' + o'qi II, and M = X(aQ/)Z(&,g/). Next we shall calculate a 
lower bound for the fidelity between \ip) and the recovered state M~ 1 Tr^ onv (|'0 / )('0 / |)(M _1 )* 
when is projected to |<tq' + 0q/) 6 Q' ®H enw after the measurement. Observe that taking 
partial trace over if env and applying M _1 to \ctqi) and Icqz + ctq,) yields the original state \ip) 
and the recovered state M~ 1 TrH cnv (|^ / )('?/' / |)(M~ 1 )*, respectively. The fidelity between \ip) 
and M~ lr Tv Hmv {\ip') (ip'\)(M~ 1 )* is not less than that between \ctqi) and \ctqi + (Tq,), because 
the fidelity does not decrease by unitary operations and taking partial trace | 12|| . We shall 
calculate a lower bound for the fidelity Fqi between |<tq/) and \o~qi + o~'q,). 



Fq' 



(<7Q'WQ')(°Q' + a Q'\ a Q' + a Q>) 



< 



{oq'\oq<){o'q> 


°'q>)- 


-K 


\ a Q'){ a Q' 


a Q') 


( a Q'\ a Q'}WQ> + 
{ a Q'\ a Q') {o'q 1 


°'q>) 


T Q'+ CX Q>) 



{°Q'\°Q'){°Q' + v'q'Wq' + °'q') 

= {<7q>W'q>) 

(cTQ' + Vq'Wq' + Vq')' 

We shall calculate a lower bound for the average of 1 — Fqi, where the average is taken 
over the measurement outcome. The following fact will be used. For a pair of vectors (a, b), 
w{a,b) denotes w(X{a)Z(b)). 
Proposition 4 

|J T Q , c{(a,6)GZ«xZ£ : w(a,b)>t}. 

Q' is an eigenspace of S 

Proof. The assertion follows from Proposition 1. ■ 
In the following calculation Q' runs through the set of eigenspaces of S. 

Y,P Q >{l-F QI ) < J2(a' QI \a' QI ) 

Q' Q 1 



E 



J2 X(c Q ,)Z(d Q ,)\<p) ® L SQfiQi \o eay ) 

(c ' 4 i)&T , 
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< E E ||x(c Q o^(4)l^)®^ Q/( r Q ,|o cnv )|| 2 

< J2 \\ X ®Z$\<P) ® ^edOir (by Proposition 4) 

«>(c,d)>t 

= E H^-OcnvH 2 . (4) 



u>(c,d)>t 

For a vector a = (ai, . . . , a n ) G Zg, let 



£(0) = ||L 0j oO E || 2 , 

£(1) = ||L ,iO £ || 2 + ||L li0 £ ;|| 2 + ULliObII 2 , 

n 

A(a) = n^)> 

i=i 

h(a) = |}{z : 7^ 0}. 



Observe that A(a) < by Assumption 3 and Lemma 2. For vectors a, b G Zg, let or (a, 
6) be the bitwise logical or of them. By these notations, for a vector a G Z2 we can see 



E ||^-0en V || 2 = A(a)<^, 

c,Jez™ 

or(c,<i)=a 

and w(c, d) = h(or(c, d)). By these observations we can rewrite Eq. (f|) as 

E ll^-Oenvf = E A(O) (5) 
«,(c,d)>i h ( s )>* 

< E p m (6) 



h(a)>t 



i=t+l \ / 



Thus 



i=t+l 



is a lower bound for the average of the fidelity between the original state and the state 
recovered by a t-error correcting quantum code of length n. 

Example 5 By Ref. || Table III] it is known that there exists a [[25, 5, 7]] code. We take it 
as an example. Then we have t = 3. At p = 0.01, the value of Eq. (|7|) is 1 — 0.000132, and 
at p = 0.001 the value of Eq. (g) is 1 - 0.127 x 10~ 7 . 



IV. CONSEQUENCES AND GENERALIZATIONS 
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A. Error-free communication is asymptotically possible 



In classical information transmission we can make the error probability arbitrary small 
by increasing the code length. The same result also holds in the quantum case. Let a be a 
real number such that 2p a < 1. Suppose that there exists a sequence of tj-error correcting 
quantum codes of length Hi such that U/rii — > a and n$ — > oo as i — > oo. The existence of 
such a sequence is guaranteed by the quantum Varshamov-Gilbert bound Theorem 2] in 
certain range of a. 

We shall consider the asymptotic behavior of Eq. (|7|): 



If t/n > a then p t+1 2 n < p(2p a ) n , which converges to as n — > oo. Thus if we use the 
sequence of quantum codes described above, we can make the average of fidelity arbitrary 
close to 1 by increasing the code length. Note that our estimate differs by factor of 2 n from 
an intuition 1 — 0(p t+1 ) of the fidelity of quantum error correction. 



The memoryless assumption is used only in Eq. We can calculate a lower bound 
for the average of the fidelity over an arbitrary channel as Eq. (||) by rewriting the unitary 
operator in a unitary representation as 



We can generalize the result to nonbinary stabilizer codes as follows. We consider q- 
ary stabilizer codes. Let H q be the g-dimensional Hilbert space and |0), \q — 1) an 
orthonormal basis of H q . Let A be a primitive g-th root of 1, for example, exp(27ri/g). We 
define a linear map C q sending |z) to \i + 1 mod q) and D\ sending \i) to X l \i) flT3| . Observe 
that Ci = <J X and D_i = o z when q — 2. 

Let T be the channel superoperator on H q , and suppose that a unitary representation of 




B. General channel 



J2 X{a)Z{b)®L^. 



C. Nonbinary codes 




We can write U as 



where Z q = {0, . . . , q — 1} and L^j is a linear operator on H E . 
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Replace the definition of p in Assumption 3 with 

p= \\ L i,jO E \\ 2 - 
(o,o)?fe(ij)ez| 

Then the lower bound Eq. (|7]) also holds for g-ary stabilizer codes. 

D. Bounded distance decoding 



In the error correction process described in Sec. |IIB| we have to find an operator M G E 
such that w(M) is minimum among operators iV G E such that NQ = Q'. The task of 
finding such M from the measurement outcome becomes computationally difficult when 
both the code length and the minimum distance are large [16]. In practice, we may give up 
finding such M if there is no operator N of weight < t' such that NQ = Q', where t' is an 



integer < t. This is a quantum analogue of the classical bounded distance decoding []I3|| . We 
shall slightly modify this bounded distance decoding and give a lower bound for the average 
of fidelity. 

Let Q, Q', \if)') and I env be as in Sec. [11 B. If there is an operator N G E such that 



NQ = Q' and w(N) < t', then let M = N. Otherwise choose an operator M G E such that 
MQ = Q' . Let the recovered state be (M _1 ® ienv)!^')- With this error correction process 
the average of the fidelity is bounded from below by 




The proof is almost the same as that of Eq. ([/]) . 

E. Nonstabilizer codes 

It seems difficult to generalize the result in this paper to nonstabilizer codes. Because in 
the error correction of nonstabilizer codes we cannot write as sum of eigenvectors of the 
measured observable as in Eq. (|3]). 

F. Entanglement fidelity 



The entanglement fidelity introduced in Refs. 0, [TTJ] should also be considered in some 
applications, and we can estimate the entanglement fidelity from the fidelity by their relation 
0, Theorem V.3]. 
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